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ABSTRACT 

In this paper, we study the filamentary structures and the galaxy alignment along 
filaments at redshift z = 0.06 in the MassiveBlack-II simulation, a state-of-the-art, 
high-resolution hydrodynamical cosmological simulation which includes stellar and 
AGN feedback in a volume of (100 Mpc//i)^. The filaments are constructed using the 
subspace constrained mean shift (SCMS; Ozertem & Erdogmus (2011) and Chen et al. 
(2015a)). First, we show that reconstructed filaments using galaxies and reconstructed 
filaments using dark matter particles are similar to each other; over 50% of the points 
on the galaxy filaments have a corresponding point on the dark matter filaments within 
distance 0.13 Mpc//i (and vice versa) and this distance is even smaller at high-density 
regions. Second, we observe the alignment of the major principal axis of a galaxy 
with respect to the orientation of its nearest filament and detect a 2.5 Mpc/h critical 
radius for filament’s influence on the alignment when the subhalo mass of this galaxy is 
between IO^Mq/H and 10^^Mq/H. Moreover, we find the alignment signal to increase 
significantly with the subhalo mass. Third, when a galaxy is close to filaments (less 
than 0.25 Mpc/h), the galaxy alignment toward the nearest galaxy group depends on 
the galaxy subhalo mass. Finally, we find that galaxies close to filaments or groups 
tend to be rounder than those away from filaments or groups. 
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1 INTRODUCTION 

Weak gravitational leasing by large-scale structure of the 
Universe, commonly known as cosmic shear, is a promis¬ 
ing technique to constrain cosmological parameters. It is 
among the key science cases of various ongoing and upcom¬ 
ing surveys, such as the Kilo-Degree Survey (KiDS^), Hyper 
Suprime Cam (HSC^), Dark Energy Survey (DES^), Large 
Synoptic Survey Telescope (LSST"*), Euclid® and the Wide- 
Field Infrared Survey Telescope (WFIRST®). Theoretically, 
it provides valuable information on both the geometry and 
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structure growth of the Universe without assumptions on the 
relationship between the luminous and dark matter; see e.g 
Hoekstra et al. (2002); Van Waerbeke et al. (2005); Jarvis 
et al. (2006); Schrabback et al. (2010); Lin et al. (2012); Hey- 
mans et al. (2012); Jee et al. (2013); Huff et al. (2014). If 
systematics can be sufficiently controlled for future missions, 
cosmic shear has the potential to be the most constraining 
cosmological probe, thus understanding and controlling sys¬ 
tematics in the weak lensing measurements needs to improve 
signihcantly over the next few years for its potential to come 
to fruition with large data volumes in these upcoming sur¬ 
veys. 

Cosmic shear is typically measured through two-point 
correlations of observed galaxy ellipticities. In the weak 
lensing regime, the observed ellipticity of a galaxy is the 
sum of its intrinsic ellipticity et and gravitational shear 
7 : Cobs ~ et -I- 7 . When the intrinsic shapes of galaxies are 
spatially correlated, these intrinsic alignment (lA) correla- 
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tions can contaminate the gravitational shear signal and lead 
to biased measurements if not properly removed or modeled. 
Since early work establishing its potential effects (Croft & 
Metzler 2000; Heavens et al. 2000; Catelan et al. 2001; Crit¬ 
tenden et al. 2001), lA has been examined throngh observa¬ 
tions (e.g., Hirata et al. 2007; Joachimi et al. 2011; Blazek 
et al. 2012; Singh & Mandelbanm 2014), analytic model¬ 
ing, and simnlations (e.g., Schneider et al. 2012; Tenneti 
et al. 2014, 2015). See Troxel & Ishak (2014); Joachimi et al. 
(2015), and references therein, for recent reviews. 

To create a fully predictive model of lA would require 
the understanding of the complex processes involved in the 
formation and evolution of galaxies and their dark matter 
halos, as well as how these processes couple to large-scale 
environment. The shapes of elliptical, pressure-supported 
galaxies are often assumed to align with the surrounding 
dark matter halos, which are themselves aligned with the 
stretching axis of the large-scale tidal field (Catelan et al. 
2001; Hirata & Seljak 2004, 2010). This tidal alignment 
model leads to shape alignments that scale linearly with fluc¬ 
tuations in the tidal field, and it is thus sometimes referred to 
as “inear alignment”, although nonlinear contributions may 
still be included (Bridle & King 2007; Blazek et al. 2011, 
2015). As for spiral galaxies, where angular momentum is 
likely the primary factor in determining galaxy orientation, 
lA modeling is typically based on tidal torquing theory, lead¬ 
ing to a quadratic dependence on tidal field fluctuations (Pen 
et al. 2000; Catelan et al. 2001; Hui & Zhang 2002; Lee & 
Pen 2008), although on sufficiently large scales, the contri¬ 
bution linear in the tidal field may become dominant. 

Filaments in the cosmic web play a key role in the 
large- scale tidal field since they are related to the gra¬ 
dient and the Hessian matrix of the matter density field 
(Hahn et al. 2007a,b; Sousbie et al. 2008; Forero-Romero 
et al. 2009; Aragon-Calvo et al. 2010; Cautun et al. 2013). 
As a result, it is expected that filaments would affect the 
shapes and alignments of their nearby galaxies. There is ev¬ 
idence for anti-alignment between the galaxy spin orienta¬ 
tion of disk galaxies and the direction of nearby filaments 
(Libeskind et al. 2013; Tempel et al. 2013; Tempel & Libe- 
skind 2013; Aragon-Calvo & Yang 2014; Dubois et al. 2014; 
Welker et al. 2014; Laigle et al. 2015; Zhang et al. 2015), 
and an alignment between galaxy shape and the orientation 
of the nearest filament (Altay et al. 2006; Zhang et al. 2009; 
Tempel & Libeskind 2013; Zhang et al. 2013). Other evi¬ 
dence also indicates that direction of a galaxy pair (Tempel 
& Tamm 2015) and the alignment of satellite galaxies to¬ 
ward their host galaxies (Tempel et al. 2015) are dependent 
on their nearby filaments. The effect of filaments on intrinsic 
galaxy shapes is worth studying since filaments are expected 
to trace the tidal fields directly and yet the current theoret¬ 
ical model for lA only includes halos but does not explicitly 
include filaments; see e.g. Schneider & Bridle (2010). 

One can, in principle, test and constrain lA models 
and marginalize over free parameters in specific models; this 
process could be improved in the quasilinear and nonlinear 
regime via the tighter priors on model parameters that we 
could obtain with a 3D filament map of the universe which 
overlaps with sources in the weak leasing fields. In this pa¬ 
per, we apply the subspace constrained mean shift (SCMS) 
filament finder to the source density in a cosmological hy¬ 
drodynamic simulation to test whether we can use filaments 


as a tool to understand the relationship of galaxy alignments 
with the large-scale density field. 

We start with a brief introduction about density ridges, 
the model for filaments we are using, in § 2 and describe 
the MassiveBlack-H Simulation we applied in § 3. Then we 
describe our main results in S 4 and our conclusions in S 5. 


2 DENSITY RIDGES AS FILAMENTS 


We detect filaments through SCMS introduced in Ozertem & 
Erdogmus (2011); Chen et al. (2015a). SCMS is a gradient- 
based method that detects filaments as “galaxy density 
ridges”. 

Essentially, filaments from SCMS are defined as follows. 
Given galaxies located at Xi, ■ ■ ■ , A„ € R®, the smoothed 
density field (also known as the kernel density estimator) is 


p{x) 



i=l 




( 1 ) 


where A is a Gaussian and b is the smoothing bandwidth 
that controls the amount of smoothing on each particle. In 
this paper, we apply 6=1 Mpc/6. Note that the density field 
p does not incorporate any boundary condition at the edge 
of the simulation box. Let g{x) — \7p{x),H{x) = Wp{x) 
be the gradient and Hessian matrix of p{x). Ridges of p{x) 
are the collection of points 


F = {x ■. Vj{x) ■ g{x) = 0, A 2 (x) < 0,j = 2,3}, (2) 

where Vj{x), \j{x) are the j-th eigenvector/value for H{x). 
See Figure 1 for an example about density ridges. More prop¬ 
erties of ridges of a function can be found in Eberly (1996); 
Ozertem & Erdogmus (2011); Genovese et al. (2014); Chen 
et al. (2014, 2015c). 

Intuitively, at a point on the ridge, the gradient is the 
same as ui(a;), the eigenvector corresponds to the largest 
eigenvalue, and the density curves downward sharply in di¬ 
rections orthogonal to the gradient. When p is smooth and 
the eigengap 


fi{x) = \i{x) - \ 2 {x) (3) 

is positive, the ridges have the properties of filaments (Chen 
et al. 2015c). That is, F decomposes into a set of smooth 
curve-like structures with high-density. Since the ridges are 
smooth, for any point x £ F, the orientation of ridges at x, 
denoted as pRidge(®), is well-defined. 


3 MASSIVEBLACK-II SIMULATION 

In this work, we detect filaments in the MassiveBlack-H 
(MBH) hydrodynamic simulation (Khandai et al. 2015), and 
study the alignments of the shapes of stellar matter compo¬ 
nents in galaxies with the filaments. MB-H is a state-of-the- 
art, high-resolution cosmological hydrodynamic simulation 
performed in a cubic periodic box of size 100 Mpc//i on a 
side. The simulations have been performed with the code 
p-GADGET, which is a hybrid version of the parallel code 
GADGET2 (Springel et al. 2005) that has been modified and 
upgraded to run on Petaffop-scale supercomputers. MBH in¬ 
cludes the physics of multiphase interstellar medium (ISM) 
model with star formation (Springel & Hernquist 2003), 
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Figure 1. Illustration of density ridges (blue curves) in a den¬ 
sity field. Each black dot is a subhalo and the big red dots are 
galaxy groups. The color from white-green-yellow-orange denotes 
the density induced by smoothing the subhalos. Note that this 
2-D figure is for illustration; our analysis is applied to 3-D data. 


black hole accretion and associated feedback processes (Di 
Matteo et al. 2008, 2012) in addition to gravity and hydro¬ 
dynamics. This simulation has been run from 2 = 159 to 
2 = 0.06 using Afpart = 2 X 1792® dark matter and gas par¬ 
ticles with a gravitational smoothing length, e = 1.85kpc/h 
in comoving units. The mass of each dark matter particle is 
moM = 1.1 X IO’^Mq/H and the initial mass of a gas parti¬ 
cle is rugas = 2.2 X lOPM q/H. The cosmological parameters 
in the simulation chosen according to WMAP7 (Komatsu 
et al. 2011) are as follows: amplitude of matter fluctuations 
as = 0.816, spectral index Us = 0.96, mass density param¬ 
eter Q,m = 0.275, cosmological constant density parameter 
Q.A = 0.725, baryon density parameter Q.b = 0.046, and 
Hubble parameter h = 0.702. 

To generate halo catalogs, the friend-of-friends (FOF) 
halo finder algorithm (Davis et al. 1985) is used with a link¬ 
ing length of 0.2 times the mean interparticle separation. 
The SUBFIND code (Springel et al. 2001) is used on the halo 
catalogs to generate subhalo catalogs. Here, subhalos are de¬ 
fined as locally overdense, self-bound particle groups which 
consist of at least 20 gravitationally bound particles. How¬ 
ever, to ensure that the measured shape of a galaxy is reli¬ 
able, we restrict our galaxy sample to those subhalos with 
at least 500 star particles. This is based on the convergence 
tests in Tenneti et al. (2014). 


3.1 Shapes of galaxies 

In this paper, we use the shapes of the stellar components 
of galaxies in MBH which are calculated in Tenneti et al. 
(2015) to quantify intrinsic alignments of galaxies. Here, we 
give a brief description of the method adopted to calculate 
these shapes. The shapes of galaxies are modeled as ellip¬ 
soids in three dimensions by the eigenvalues and eigenvectors 



Figure 2. Example of /rc', fip, dc, dp. The big red dot denotes a 
galaxy group, the blue curve is the filament, and the small black 
points are galaxies. For a given point (galaxy), the brown arrow is 
fiQ and the purple arrow is /ip. The length of the dotted brown 
line is dc and of the dashed purple line is dp. Note that the 
projected point onto the filament is the pink point. 


obtained from the eigen-decomposition of the reduced iner¬ 
tia tensor. We consider all the star particles belonging to the 
subhalo in the calculation of the inertia tensor. The eigenvec¬ 
tors of the inertia tensor are denoted as /ri,/t 2 , and /ta with 
the corresponding eigenvalues being Xa,\b, and Ac, where 
Aa > Ai, > Ac. The lengths of the principal axes of the ellip¬ 
soid (a, b, c), which are represented by these eigenvectors, are 
given by square roots of the eigenvalues (vX(, \/A6, vX)). 
We can now define the intermediate-to-major axis ratio (q) 
and the minor-to-major axis ratio (s) as 


Note that an iterative procedure is adopted in calculat¬ 
ing the shapes such that the axis ratios converge. For ad¬ 
ditional details regarding the iterative method, see Tenneti 
et al. (2015). 

For each galaxy, we have two intrinsic quantities: the 
number of star particles n and the total mass of this sub¬ 
halo M (subhalo mass). Throughout this paper, we use only 
subhalo mass. 


3.2 Large-Scale Structure 

In our analysis, we define galaxy groups as those subha¬ 
los with subhalo mass > This gives 315 galaxy 

groups, which corresponds to a comoving number density 
3.15 X 10~'*(h/Mpc)®. Let C and F denote the collection of 
all groups (of galaxies) and filaments. Given any galaxy lo¬ 
cated at X {x is the 3D position in the box), let Ilc{x) € C 
and nF(x) & F he the ‘projected’ point onto the nearest 
group/filament. Namely, nc(2:) (or nF( 2 ;) respectively) is 
the point in C (or F) whose distance to x is smallest among 
all points in C (or F). In general, this projection is unique 
(although some regions may have non-unique projection, the 
probability that a galaxy falls within these regions is 0). We 
then define the associated direction to the nearest group and 
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Notation 

Definition 

Type 

Remark 

n 

Number of star particles 

Integer 

> 500 

M 

Halo mass 

Scalar 


dc 

Distance to the nearest galaxy group 

Scalar 


dp 

Distance to the nearest filament 

Scalar 


Me 

Direction to the nearest galaxy group 

Vector 


lip 

Orientation of the nearest filament 

Vector 


Q 

First and second principal axis ratio 

Scalar 

6 [0,1] 

s 

First and third principal axis ratio 

Scalar 

6 [0,1] 

Ml 

First principal axis 

Vector 


M2 

Second principal axis 

Vector 


M3 

Third principal axis 

Vector 



Table 1. Table of quantities associated with a galaxy. 


to the nearest filament as 


nc(a:) — X 

|nc(®) - x\ 

(5) 

Ridge (np (x) ) . 

(6) 


Note that is the direction toward the nearest galaxy 
group while /tf is the direction of the nearest filaments. 
We also define the projection distances dc{x) = |nc'(a;) — x\ 
and dp(x) = |nF(ai) — 2;|. Namely, for a galaxy at position x, 
Hc(x) is the direction from this galaxy to the nearest galaxy 
group. And ^,f(x) is the orientation of the nearby filament. 
Figure 2 illustrates all the above quantities in a 2-D picture. 


4 RESULTS 

4.1 Agreement of Filaments from Galaxies versus 
Dark Matter Particles 

To establish the promise of recovering filaments tracing dark 
matter particles using only galaxies, we first study the sim¬ 
ilarity for filaments obtained by using galaxies only versus 
using dark matter particles. The dark matter particles are 
from a random subsample with size 400,000, which is ap¬ 
proximately the same order as the number of subhalos we 
have (~ 500,000). In Figure 3, we show the difference be¬ 
tween the two filament catalogues. Since both catalogues 
contain points on the filaments, the difference is computed 
using projected distance between points from two catalogues 
(Chen et al. 2015b). To be more specific, let To and Fg be 
the collection of points on the filaments from dark matter 
particles and from galaxies. Then for each point in To, we 
find its nearest distance to any member of Fa- This gives 
a distance value (also known as projected distance) to each 
point in Fo. Similarly, for every point in Ta, we can find 
the corresponding projected distance to J-d so that we will 
assign a distance value to each element of J-q- The distri¬ 
bution of all the distance values in and J-q provides 
the information about similarity between these two filament 
catalogues. If they are similar, this distribution will be con¬ 
centrated around 0. On the other hand, if the two catalogues 
are different, this distribution will spread out. Note that we 
have tried a random sample with size 800,000 and the result 
remains similar. 

According to Figure 3, we see that the two catalogues 
are similar to each other; more than 50% of the points on 


Difference of filaments using gaiaxies versus DM 



Figure 3. Agreement of filaments detected using galaxies ver¬ 
sus using dark matter particles. We apply our filament finder to 
both galaxies and dark matter particles and then compare the dif¬ 
ference between the two corresponding filament catalogues. This 
picture shows the distribution of difference in filament position. 
The four color bars denote the quantile values of 25%, 50%, 75% 
and 90%. The difference for filaments recovered from galaxies and 
dark matter particles is small. 


filaments differ at a distance less than 0.13 Mpc//i and about 
90% of them differ at the distance less than 0.5 Mpc//i. 

We also study the agreement for the two catalogues at 
different densities. For each point on the filaments, we assign 
a density value to them by using the distance to the 1000-th 
nearest dark matter particles (Figure 4). The 1000 is cho¬ 
sen arbitrarily; we have tried other values and the result 
remains similar. The inverse of this distance is an effective 
density (at overdense regions, this distance is small while at 
void region, this distance is large). We then rank points on 
filaments according to their densities and plot the distribu¬ 
tion of difference under different density ranking categories. 
The high-density filaments are those whose density ranking 
is among the top 10% (red) and top 25% (orange). The low 
density filaments are those whose density ranking is in the 
bottom 25% (green) and bottom 10% (blue). The result is 
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Density: dark matter particles 



Figure 4. Agreement of filaments detected using galaxies versus using dark matter particles under different environmental densities. 
We rank points on filaments according to the densities of dark matter particles (using densities from galaxies yields a similar result). 
Then we plot distributions of filament differences under different density regimes. Filaments at high-density regions are colored by red 
and orange (density within the top 10% and 25%) whereas low density filaments are colored by green and blue (density in the bottom 
25% and 10%). The black curve denotes the distribution using all galaxies (which is the distribution of Figure 3, note that we use the 
different bin size). We see a clear pattern that the two filament catalogues are similar at high-density regions and get less similar when 
the density decreases. 


given in Figure 4. We observe that the distribution of differ¬ 
ence for both catalogue is small, and the difference is even 
smaller for Hlaments reside in high-density regions. The rea¬ 
son why the difference is density-dependent is because galax¬ 
ies are a biased tracer for the matter density held so that 
galaxies trace the large-scale structure structures better at 
high-density regions. Note that the result remains similar 
when we replace dark matter density by galaxy density. 

To conclude, Figure 3 and Figure 4 show that the hl¬ 
aments from galaxies and the hlaments from dark matter 
particles are similar and this similarity is stronger for high- 
density regions. For the following analysis on the galaxy 
alignment, we use the hlaments from galaxies. 


4.2 Principal Axes Alignment to Filaments 

We hrst examine the correlation between principal axes and 
Hf, the orientation of the nearest hlament of a galaxy. Note 
that since both and —fj,j are describing the same axis, we 
take the absolute value of inner product between fij and hf- 
We compute the alignment signal • /tfI for galaxies 
close to hlaments < 1 Mpc/h) but not close to groups 
{dc > 2 Mpc//i). Note that • ^f\ ~ 1 means that the j-th 
principal axis align along /if while l/ij ■ ^f| ~ 0 indicates 
that they are misaligned. In Figure 5, we show the excess 
probability compared to the distribution of inner product 
between two random unit axes. We hnd that the alignment 
signal for hrst principal axis (blue curve) is very signihcant 
(with KS-test p-value being 4 x 10“^®). By comparing the 
blue curve in both the left and right side of Figure 5, we 
see that the distribution of the alignment is concentrated 
at values « 1. This means we have more aligned galaxies 


compared to two galaxies with random orientations. The 
trend is reversed for the minor principal axis (with KS-test 
p-value 1 X 10“^^). Thus, we conclude that the major (hrst) 
principal axis for a galaxy tend to align along hlaments while 
the minor axis tends to misalign to the hlaments. 

To study how the subhalo mass of a galaxy affects the 
alignment signal, we further partition galaxies by their mass 
into three mass-groups 

10®Mo//i < M < W^°MQ/h, 

< M < 1O“M0//i, (7) 

lO“M0//i < M < W^^MQ/h 

and plot the average alignment signal • fiF\) within each 
mass bin. Note that the above three mass-groups contain 
88 % of all subhalos we have (about 11 % of the subhalos have 
subhalo mass M < IO^Mq/H and 0.6% have subhalo mass 
M > lO^^M 0 /h). The result is given in the hrst two panels 
in Figure 6 . The left (hrst) panel shows a clear patten that 
the average alignment signal for major axis increases as the 
subhalo mass increases (this slope has a du signihcance) and 
the effect of misalignment for minor axis is also augmented 
for massive galaxies. Altay et al. (2006); Zhang et al. (2009); 
Tempel & Libeskind (2013); Zhang et al. (2013) also report 
a similar mass dependency on the alignment signal. 

The middle (second) panel of Figure 6 shows the same 
analysis as the left panel but at a different distance range 
to hlaments. We now focus on galaxies whose distance to 
hlaments is within 1 Mpc//i < dp < 2 Mpc/h. That is, the 
middle panel shows the alignment signal for galaxies that 
are mildly close to hlaments (but not too close). Our result 
shows that, unlike galaxies being very close to hlaments (left 
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Figure 6. Alignment signals \{{ij • fiF)\ versus subhalo mass and djr, the distance to filaments. Note that we remove galaxies that are 
close to groups (i.e. all galaxies satisfy dc > 2 Mpc/h). Left: We show how alignment signals change as a function of subhalo mass at 
dF < 1 Mpc/h. The slope of \{fii • fiF)\ versus djr is indicated in the panel. We observe a significant (4 (t) linear effect from subhalo 
mass on the alignment of first principal axis along the nearby filament orientation. Middle: The same analysis as left panel but with 
different range for distance to filaments (1 Mpc/h < dF < 2 Mpc/h). Compared with left panel, we still see an alignment signal but 
the mass-dependency disappears. Right: The alignment signal for the first principal axis, (|/ii • a function of dF- Note that we 

display the significance on top of the alignment value for each mass bin. The alignment signal is significant when dF < 2.5 Mpc/h (left 
side of the orange line). 


dp < 1 Mpc/h, dc > 2Mpc/h 



Figure 5. Distribution of alignment signals \{iij ■ /if)I - Note that 
we only consider galaxies away from groups {dc > 2 Mpc/h) 
but close to filaments {dp < 1 Mpc/h). We observe a significant 
increase for |(/ii ■ /tf}! ~ 1- This shows that the major principal 
axis tend to align along the orientation of the nearby filament. 
The 5 in the Y-axis is the excess probability compared to the 
distribution of random angle. The error bars are calculated from 
the bootstrap. The number below shows the sample size for each 
bin. 


panel), the alignment signal is not mass-dependent (only 
0.2cr). 

On the right panel of Figure 6, we plot the average align¬ 
ment signal for the first principal axis at different distance 
bins to show how the distance affects the alignment signal. 
Note that here we only focus on galaxy with subhalo mass 
between [IO^Mq/Zi, IO^^Mq/Zi]. As expected, when the dis¬ 


tance to filaments is small (less than 2.5 Mpc/h), the align¬ 
ment signal is significant (> 2a). However, when the dis¬ 
tance is above 2.5 Mpc/h, we do not observe any significant 
alignment signal. This result is consistent with the left and 
middle panels of Figure 6 since the two panels both satisfy 
dp < 2.5 Mpc/h. 

Finally, we focus on the alignment signal for first prin¬ 
cipal axis and study how the shape {q and s) influences the 
signal. We partition galaxies according to their q and s. The 
cut on q is 0.83 and the cut on s is 0.68; both cuts are chosen 
by the median value for all galaxies. The result is given in 
Figure 7. The left panel shows the alignment signal for the 
major axis under the cut on q. We find that the two groups 
do not have significant difference (significance is only 0.42 ct). 
The right panel shows the same analysis under the cut on 
s. We see a significant difference (3.36 ct) for the two groups 
of galaxies. Namely, an elongated galaxy (with smaller s) 
tends to align more along nearby filaments compared to a 
spherical one. This is not an effect due to measurement error 
since in the simulations, all the galaxies have > 500 particles 
and the shapes and orientations are well-measured even in 
the sample that are closer to round. 

4.3 Principal Axes Alignment to Groups 

We also study how filaments influence the galaxy alignment 
toward the nearest galaxy groups. We focus on the inner 
product between the major (first) principal axis and the di¬ 
rection to the nearest galaxy groups. Namely, we test how 
the alignment signal ■ /ic| is different for galaxies within 
a filament versus outside filaments. To classify a galaxy as 
being inside a filament or out of filaments, we call a galaxy 
in filaments if their distance to filament dp is less than a 
distance threshold r. Otherwise, we say this galaxy is out 
of filaments. Note that we only consider galaxies that are at 
least 1.00 Mpc/h away from the nearest group to make sure 
this galaxy is not in this group. 
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Alignment of the First Principal Axis Alignment of the First Principal Axis 




Figure 7. The correlation between shapes and the alignment signal. Note that we focus on the alignment signal of first principal axis. 
We study how a galaxy’s shape is correlated with \{fii • We split the galaxies into groups according to their shapes, measured by 

the q and s. The separation cuts on q and s are chosen by the median of q and s for galaxies with dp < 1 Mpc//i and dc > 2 Mpc//i . 
Left: We separate galaxies into two groups by q\ more elliptical galaxies [q < 0.83, orange) and more spherical galaxies (q > 0.83, red). 
There is no significant difference between these two groups (0.42(t). Right: We separate galaxies into two groups by s: more elliptical 
galaxies (s < 0.68, green) and more spherical galaxies (s > 0.68, blue). In this case, we observe a significant difference between the two 
groups (3.36cr). An elliptical galaxy seems to align more along the filaments compared with a spherical one. 
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Figure 8. The effect for filaments on |(/ri ■ /rc}|, the galaxy alignment along ‘the direction toward the nearest group’. We focus on the 
galaxies with dc > 1 Mpc//i to make sure they are not inside a group. We then partition galaxies into two types: in filaments and out 
filaments by imposing a threshold r on the distance to filaments (see § 4.3 for more details). Left: t = 0.25 Mpc//i. Namely, a galaxy 
whose distance to filaments is less than 0.25 Mpc/h will be considered within filaments. Otherwise this galaxy is classified as out of 
filaments. Middle: r = 0.50 Mpc/h. Right: r = 1.00 Mpc/h. The red color denotes |(/ri ■ /re) I for galaxies that are inside a filament 
and the blue color is the same but for galaxies out of filaments. Note that we also attach the significance for the slope using the linear 
regression. 


The result is given in Figure 8. We consider three dif¬ 
ferent distance thresholds: 0.25 Mpe/h (left), 0.50 Mpe/h 
(middle), and 1 Mpe/h (right). In every panel, we observe 
that, for galaxies in filaments (red curves), the alignment sig¬ 
nal significantly depends (the significance ranges from 2.4 ~ 
3.Oct) on the subhalo mass for the galaxy. On the other hand, 
the alignments signal for galaxies outside filaments (blue 
curves) are significant only when we set r = 0.25 Mpe/h 
but this dependency disappears for other cases. Moreover, 
when we increase the threshold r, the slope (the linear effect 
from mass on galaxy alignments toward groups) for galaxies 
within filaments decreases but the slope for galaxies out of 
filaments increases. This suggests that the majority of the 
mass-dependency is from galaxies with dp < 0.25 Mpc//i. 


4.4 Shape Analysis 

Now we study the correlation between the shape of galaxies 
and their proximity to galaxy groups and filaments. In par¬ 
ticular, we focus on q for each galaxy since it measures the 
ratio to the length of first and second principal axis and is 
a useful summary statistic for the ellipticity. 

We partition galaxies according to their subhalo mass 
into four groups: 

10^° °Mo//i < M < lO^° ®M0//i, 
lO^° ®M0//i < M < 1O“ “M0 //i, 
lO^^ °M0//i < M < lo“-^M0//i, 

1O“ '‘M0//i < M < 1O^^-“M0//i. 

Note that this partitioning is different from equation (7) 
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dc < 5 Mpc/h 


dp < 5 Mpc/h, dc > 2 Mpc/h 



Figure 9. Impact from subhalo mass on the shapes of a galaxy. Left: the result for galaxies close to groups (dc < 5 Mpc/h). Right: 
the result for galaxies close to filaments but not close to groups (dp < 5 Mpc/h and dc > 2 Mpc/h). The two distance ranges will be 
used in Figure 10. In both cases, we observe a strong impact on both q, s from subhalo mass. 


because here we separate galaxies by subhalo mass to control 
the effect from subhalo mass on q. Therefore, we need a 
fine mass bin to reduce the mass-dependency on q. Figure 9 
shows how (q) and (s) change as the subhalo mass changes. 
This pattern corroborates what is observed in Tenneti et al. 
(2014). Note that Schneider et al. (2012) observed a similar 
result using dark matter halos and dark matter halo shapes. 

Figure 10 displays how (q) changes as dp or dc changes. 
We stack galaxies around filaments (and groups) and then 
partition galaxies according to their distance to filaments 
(and groups) into 10 bins with a bin width of 0.5 Mpc/h (so 
we only consider distance to filaments/groups less than 5 
Mpc/h). For galaxies stacking around filaments, we have an 
additional constraint dc > 2 Mpc/h to eliminate the effect 
from galaxy groups. For each bin, we also show the average 
subhalo mass. 

Overall, we find two cases that lead to a change in the 
(q) after adjusting the effect from subhalo mass. The first 
case can be found at the top right and bottom left panel 
of Figure 10. We observe a decreasing pattern as dp or dc 
increases. In other words, for galaxies close to filaments or 
groups, we expect their q to be higher (more rounded galax¬ 
ies). In both panels, the subhalo mass within each bin does 
not change too much so that this effect cannot be explained 
by the mass difference. The second case is on the difference 
between groups and filaments. In every panel of Figure 10, 
we find that, for galaxies at the same distance to either fil¬ 
aments or groups, (q) is always higher for galaxies close to 
groups. 

An explanation to the above cases is due to the differ¬ 
ence in the ratio of elliptical galaxies. If elliptical galaxies 
are more likely to be in/near groups and filaments and if 
they are usually rounder than disk galaxies, then we would 
expect to see that galaxies close to groups and filaments are 
on average rounder than those that are far away. And if this 
effect is stronger for groups (perhaps due to a stronger tidal 
field), then we should see a higher q for galaxies that are 
closer to groups. 

Note that in the upper left panel (IO^^'^Mq/Zi < M < 
Mq/K), we find a decreasing trend for the average sub¬ 
halo mass when the distance to large-scale structure is in¬ 
creasing. Moreover, at the same distance, galaxies close to 


groups tend to be lower mass compared with those close 
to filaments. Namely, one will find more low-mass galaxies 
(with subhalo mass less than IO^^'^Mq/Zi) at strong tidal 
force regions. We do not find a similar patten in other mass 
range, so this phenomena occurs only for light galaxies. A 
possible reason is that, in a region with strong tidal ef¬ 
fect, matter concentrates around the large-scale structures 
so that there will be many newly formed galaxies. These 
young galaxies are generally light, which decreases the aver¬ 
age mass. 


5 CONCLUSION 

In this paper, we study the correlation between galaxy ori¬ 
entation and filaments by applying the density ridge model 
to the smoothed particle hydrodynamic simulation. A short 
summary for what we have found is as follows: 

(i) We demonstrate that filaments constructed using density 
ridges from galaxy density field are similar to filaments con¬ 
structed from dark matter density field. This analysis, along 
with the fact that density ridges are statistically consistent 
in both theory (Chen et al. 2015c) and simulation (Chen 
et al. 2015a), suggests that we might be able to reconstruct 
filaments of dark matter field by using only galaxies. How¬ 
ever, this simulation goes down to quite low subhalo mass 
compared to most observations. More analysis needs to be 
done in future for real observations. 

(ii) We find that galaxies align along with their nearby fila¬ 
ments and this alignment is stronger for massive galaxies 
when the galaxy is close to filaments (dp < 1 Mpc/Zi). This 
corroborates the past literature (Altay et al. 2006; Zhang 
et al. 2009; Tempel & Libeskind 2013; Zhang et al. 2013). 
Moreover, we discover that the alignment signal has a sphere 
of influence of 2.5 Mpc/Zi in radius. Namely, the principal 
axes of a galaxy whose distance to the nearest filaments is 
less than 2.5 Mpc/h would be influenced by the nearby fila¬ 
ments. This critical radius has not been previously reported 
in the literature. 

(iii) We also observe an interesting pattern: the galaxy align¬ 
ment toward the nearest galaxy group is mass-dependent 
when this galaxy is near a filament within 0.25 Mpc/Zi. For 
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Figure 10. Influence from groups and filaments on the shape of galaxies. We display how (q) changes as the distance to filaments (blue) 
and groups (brown) changes. Note that we partition galaxies according to their subhalo mass to adjust the mass-dependency on (q). The 
number near (q) at each bin is the average subhalo mass within that bin. 


galaxies close to a filament, massive galaxies tend to align 
along the direction toward the nearest group. For galaxies 
away from filaments, subhalo mass has little impact on the 
alignment toward a group. 

(iv) Finally, we find that the ellipticity of a galaxy is correlated 
with its distance to groups and filaments when the subhalo 
mass is between and Mq. A galaxy that is 

close to either a group or a filament tends to be less elliptical 
(more rounded). 

The galaxy-filament alignment observed in this paper 
can be explained by the intrinsic alignment model and the 
relationship between tidal field and filaments. The intrin¬ 
sic alignment model links the galaxy principal axes to the 
tidal field (Catelan et al. 2001; Hirata & Seljak 2004, 2010) 


and the tidal field is often explicitly or implicitly used to 
define a filament (Hahn et al. 2007a,b; Forero-Romero et al. 
2009; Cautun et al. 2013). Therefore, the principal axes and 
the shape of a galaxy are expected to be correlated with its 
nearby filament. Based on the galaxy-filament interaction we 
have found, we expect the future survey on intrinsic align¬ 
ment and weak leasing to be improved if one incorporates 
the effect from filaments. This might be a way to extend the 
current linear or non-linear alignment model (Bridle & King 
2007; Blazek et al. 2011, 2015) since the alignment along 
filaments can be viewed as a new type of non-linear effect. 
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